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COMPUTER  PROGRAM 

NARFET  -  WIND-WAVE  GENERATION  ON  RESTRICTED  FETCHES 

PURPOSE :  The  purpose  of  this  technical  note  is  to  predict 

significant  wave  height,  peak  period,  and  mean  direction  for 
narrow  or  restricted  fetches. 

BACKGROUND :  Wind-wave  generation  in  lakes,  rivers,  bays,  and 

reservoirs  is  generally  limited  by  the  geometry  of  the  water 
body,  which  is  often  very  irregular.  Most  approaches  to  this 
problem  consider  wave  generation  only  in  the  direction  of  the 
wind  with  fetch  lengths  averaged  over  small  arcs  (Shore 
Protection  Manual  (SPM)  1984)  or  large  arcs  (Saville  1954). 
Donelan  (1980)  proposed  wave  generation  on  fetch  lengths  in  off- 
wind  directions  with  reduced  wind  forcing  (reduced  by  the  cosine 
of  the  angle  between  the  off-wind  and  wind  directions)  for  the 
Great  Lakes.  The  NARFET  model  (Smith  1989)  is  based  on  the 
Donelan  concept,  allowing  wave  generation  in  off-wind  directions. 
Smith  developed  expressions  for  wave  height  and  period  as  a 
function  of  fetch  geometry  and  wind  speed  based  on  I  Inear 
regressions  of  wave  data  collected  on  Puget  Sound  (Washington), 
Fort  Peck  Reservoir  (Montana),  Denison  Reservoir  (Texas),  and 
Lake  Ontario.  The  mean  wave  direction  is  determined  by 
maximizing  the  wave  period.  These  equations  differ  from  those 
given  by  Donelan  which  were  developed  for  the  longer,  more 
regular-shaped  fetches  of  the  Great  Lakes.  The  NARFET  model  is 
quick  and  inexpensive  (runs  on  a  PC),  yet  considers  the 
complexity  of  fetch  geometry. 

ASSUMPTIONS: 

1.  Waves  are  locally  generated  and  f etch- I im i ted . 

2.  Water  depths  across  the  fetch  are  deep  (depth  is  greater 
than  half  the  wave  length. 

3.  Wind  speed  and  direction  are  steady  (spatially  and 
tempora lly). 

I NPUT :  The  input  to  NARFET  describes  the  fetch  geometry  and  the 
wind  forcing.  The  program  accepts  Interactive  responses  to 
questions.  Responses  are  either  numeric  or  alphabetic. 

Alphabetic  responses  are  shortened  to  one  letter  abbreviations 
shown  In  parentheses.  Capital  letters  must  be  used.  When  a  file 
name  is  requested,  the  number  of  characters  in  the  name  Is 
limited  to  eight.  The  input  required  includes: 

1.  Fetch  Geometry.  The  fetch  geometry  is  described  by 
radial  fetch  lengths  measured  from  the  shoreline  to  the  point  of 

interest  at  even  ana le  increments.  The  user  specify  the 
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angle  (measured  clockwise  from  north)  of  the  first  radial  length, 
the  total  number  of  radials  (for  some  locations  it  is  not 
necessary  to  input  a  full  360  degree  arc  -  unspecified  lengths 
are  set  to  zero),  and  the  angle  Increment  between  radials.  Then, 
individual  radial  lengths  are  entered  Interactively.  The  program 
al lows  the  user  to  review  radial  length  entries  and  correct 
errors.  NARFET  Internally  interpolates  fetch  lengths  at  1  degree 
increments  around  the  entire  360  degree  arc.  Then  the  program 
averages  fetch  lengths  over  15  degree  arcs  centered  on  each  one- 
degree  increment.  These  fetch  lengths  are  used  to  calculate  wave 
conditions.  The  option  is  given  to  write  this  information  to  a 
file  for  future  computer  runs.  A  file  name  is  requested  when 
this  option  is  chosen. 

2.  Wind  Forcing.  Wind  forcing  is  represented  by  wind 
speed,  direction,  and  duration  over  the  water  body.  Wind  fields 
are  distorted  by  friction  effects,  so  the  measurement  elevation, 
the  boundary  layer  stability,  and  the  measurement  location  (over 
land  or  over  water)  are  also  needed  to  adjust  the  wind  speed  to 
standard  conditions.  The  simplified  corrections  to  wind  speed 
used  in  NARFET  are  based  on  these  three  factors.  The  correction 
methods  are  given  in  the  SPM.  Wind  speeds  are  adjusted  to  the 
standard  10-m  elevation.  The  air-sea  temperature  difference 
represents  the  boundary  layer  stability.  If  the  temperature 
difference  Is  unknown,  the  SPM  recommends  a  correction  factor  of 
1.1  (unstable  condition).  This  correction  factor  Is  equivalent 
to  an  air-sea  temperature  difference  of  approximately  -3  degree 
C.  Overland  wind  conditions  differ  from  overwater  conditions 
because  of  Increases  in  surface  roughness  over  land.  An 
additional  correction  is  made  if  winds  are  based  on  overland 
measurements.  An  option  is  also  available  to  correct  winds  for 
nonconstant  coefficient  of  drag  (SPM).  This  correction  was  not 
used  in  developing  the  model,  but  current  Corps  of  Engineers 
guidance  recommends  the  correction.  The  duration  input  is  used 
to  check  if  wave  generation  is  limited  by  duration.  The  program 
does  not  convert  very  short  duration  wind  observations  (e.g. 
fastest  mile  wind  speeds)  to  longer  durations. 

PROGRAM  output :  The  primary  output  of  NARFET  is  the  significant 
wave  height,  peak  period,  and  mean  direction  cor  respond i ng  to  the 
input  wind  forcing.  The  program  also  recaps  the  input  wind 
conditions  and  states  whether  the  solution  is  f etch- I im i ted , 
durat i on- I  I m i ted ,  or  fully-developed. 

SAMPLE  PROBLEM :  For  central  Puget  Sound  (Figure  1),  find  the 
significant  wave  height,  peak  wave  period,  and  mean  wave 
direction. 

Given:  wind  speed  =  15  m/s  at  10-m  elevation  over  water 

wind  direction  =  200  degrees 
wind  duration  =  5  hours 
air-sea  temperature  difference  unknown 
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Figure  1.  Southern  Puget  Sound  Geometry. 

Underlined  values  are  user's  input.  Note  the  use  of  single¬ 
letter  responses. 

NARFET 

****************************** 
PROGRAM  NARFET 


THIS  PROGRAM  CALCULATES  DEEPWATER  WAVES  FOR  RESTRICTED  FETCHES 
BASED  ON  WIND  SPEED,  WIND  DIRECTION,  AND  FETCH  GEOMETRY 
****************************** 

DO  YOU  WISH  TO  ENTER  FETCH  GEOMETRY  ( I ) NTERACT I VELY 
OR  FROM  A  (F) I LE? 


FETCH  GEOMETRY  IS  DETERMINED  BY  I NPUT I NG  RADIAL 
LENGTHS  MEASURED  FROM  THE  POINT  WHERE  YOU  WANT  WAVE 
INFORMATION  TO  THE  LAND  BOUNDARY  OF  THE  WATER  BODY 

INPUT  THE  ANGLE  INCREMENT  BETWEEN  RADIAL  MEASUREMENTS  (DEG) 

INPUT  THE  DIRECTION  OF  THE  FIRST  RADIAL  WITH  RESPECT  TO 
THE  LOCATION  OF  INTEREST  (IN  DEGREES  MEASURED  CLOCKWISE 
FROM  NORTH) 

1-3B., 

INPUT  THE  NUMBER  OF  RAD  I ALS 

1A 
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INPUT  UNITS  OF  RADIAL  LENGTHS:  ( K ) I LOMETERS 

(N)AUTICAL  MILES 
K 


(F)EET,  (M) I LES ,  OR 


1  NPUT 
10.24 

RADIAL 

LENGTH 

FOR 

138.0 

DEG 

1  NPUT 
9.75 

RADIAL 

LENGTH 

FOR 

144.0 

DEG 

INPUT 
7  ■  8.8 

RADIAL 

LENGTH 

FOR 

150.0 

DEG 

1  NPUT 

8  t  2J- 

RADIAL 

LENGTH 

FOR 

156.0 

DEG 

1  NPUT 
8,62 

RADIAL 

LENGTH 

FOR 

162.0 

DEG 

1  NPUT 
35.68 

RAD  1 AL 

LENGTH 

FOR 

168.0 

DEG 

INPUT 

28.24 

RADIAL 

LENGTH 

FOR 

174.0 

DEG 

1  NPUT 

2Q..2Q 

RADIAL 

LENGTH 

FOR 

180.0 

DEG 

INPUT 

16.01 

RADIAL 

LENGTH 

FOR 

186.0 

DEG 

1  NPUT 
12.2,6 

RAD  1 AL 

LENGTH 

FOR 

192.0 

DEG 

INPUT 
8_.  21 

RADIAL 

LENGTH 

FOR 

198.0 

DEG 

1  NPUT 
8,U 

RADIAL 

LENGTH 

FOR 

204.0 

DEG 

1  NPUT 
7.48 

RADIAL 

LENGTH 

FOR 

210.0 

DEG 

INPUT 

5.20 

RADIAL 

LENGTH 

FOR 

216.0 

DEG 

1  NPUT 
5.0  8 

RAD  1 AL 

LENGTH 

FOR 

222 . 0 

DEG 
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RECAP 

OF 

1  NPUT 

ANGLES 

AND  RAD 

1  AL 

LENGTHS 

ANGLE 

= 

138.0 

RADIAL 

LENGTH 

— 

10.24 

ANGLE 

= 

144.0 

RADIAL 

LENGTH 

= 

9.75 

ANGLE 

= 

150.0 

RADIAL 

LENGTH 

= 

7.88 

ANGLE 

= 

156.0 

RAD  1 AL 

LENGTH 

= 

8.21 

ANGLE 

= 

162.0 

RADIAL 

LENGTH 

= 

8 . 62 

ANGLE 

= 

168.0 

RAD  1 AL 

LENGTH 

= 

35.68 

ANGLE 

= 

174.0 

RAD  1 AL 

LENGTH 

= 

28 . 24 

ANGLE 

= 

180.0 

RADIAL 

LENGTH 

= 

20.20 

ANGLE 

= 

186.0 

RAD  1 AL 

LENGTH 

= 

16.01 

ANGLE 

= 

192.0 

RAD  1 AL 

LENGTH 

= 

12.76 

ANGLE 

= 

198.0 

RADIAL 

LENGTH 

= 

8.21 

ANGLE 

= 

204 . 0 

RAD  1 AL 

LENGTH 

= 

8.17 

ANGLE 

= 

210.0 

RAD  1 AL 

LENGTH 

= 

7.48 

ANGLE 

= 

216.0 

RADIAL 

LENGTH 

= 

5.20 

ANGLE 

= 

222 . 0 

RADIAL 

LENGTH 

= 

5.08 

HOW  MANY  VALUES  DO  YOU  WISH  TO  CHANGE? 

(ENTER  0  FOR  NONE) 

Q. 

DO  YOU  WISH  TO  SAVE  FETCH  GEOMETRY  FOR 
FUTURE  RUNS?  (Y  OR  N) 

Y 

ENTER  FILE  NAME  (MAX  OF  8  CHARACTERS)  TO 
SAVE  FETCH  GEOMETRY 
PUGS ■ DAT 

INPUT  UNITS  OF  WIND  MEASUREMENT  ELEVATION:  (M)ETERS  OR  (F)EET 

M. 

INPUT  WIND  MEASUREMENT  ELEVATION 

JJLu 

IS  THE  OBSERVATION  LOCATION  OVER  WATER  (W)  OR  LAND  (L)? 

)tL 

INPUT  UNITS  OF  AIR-SEA  TEMPERATURE  DIFFERENCE:  DEGREES  (C)  OR 

Q 

INPUT  UNITS  OF  WIND  SPEED:  (M) ETERS/SEC ,  ( F ) EET /SEC ,  (K)NOTS, 

OR  MILES/HOUR  (N) 

INPUT  AIR  -  SEA  TEMPERATURE  DIFFERENCE 


INPUT  WIND  SPEED,  WIND  DIRECTION  (DEG),  AND  DURATION  (HR) 
15. .200. .5, 
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INPUT  CONDITIONS: 

ADJUSTED  WIND  SPEED  (M/S)  =  22.7  (  22.7  INPUT  UNITS) 

WIND  DIRECTION  (DEG)  -  200.0 

DURATION  ( HR ) =  5.0 

AIR-SEA  TEMP  DIF  (DEG  C)  =  -3.0 


SIGNIFICANT  WAVE  HEIGHT  (M) 
SIGNIFICANT  WAVE  HEIGHT  (FT) 
PEAK  WAVE  PERIOD  (S)  =  4.7 

MEAN  WAVE  DIRECTION  (DEG)  = 


173 


1 .6 
5.3 

0 


DURATION  LIMIT  (HR)  =  2.9 


FETCH  LIMITED  CONDITIONS 

DO  YOU  WANT  TO  RUN  ANOTHER  WIND  CONDITION? 

H 

RUN  COMPLETE 
FORTRAN  STOP 


AVA I  LAB ILITY:  A  FORTRAN  listing  of  NARFET  is  available  from  the 
Coastal  Engineering  Research  Center,  Coastal  Oceanography  Branch. 
The  program  Is  internally  documented  via  comment  cards,  and 
detailed  information  about  the  development  of  the  model 
by  Smith  (1989).  For  further  information,  contact 
Smith  at  (601)  634-2079. 


is  given 
Jane  McKee 


REFERENCES: 

Done  I  an,  M . A .  1980.  "Similarity  Theory  Applied  to  the 

Forecasting  of  Wave  Heights,  Periods,  and  Directions," 

Proceedings  of  the  Canadian  Coastal  Conference,  National  Research 
Council,  Canada,  pp.  47-61. 

Saville,  T.  1954.  "The  Effect  of  Fetch  Width  on  Wave 
Generation,"  Technical  Memorandum  No.  70,  Beach  Erosion  Board, 
Washington,  DC. 


Shore  Protection  Manual.  1984.  US  Army  Engineer  Waterways 
Experiment  Station,  Coastal  Engineering  Research  Center,  US 
Government  Printing  Office,  Washington,  DC. 


Smith,  J.M.  1989(in  preparation).  "W 
Restricted  Fetches,"  US  Army  Waterways 
Coastal  Engineering  Research  Center,  V 


ind-Wave  Generat 
Experiment  Stat 
Icksburg,  MS. 


I  on  on 
i  on , 


6 


